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Monte Carlo simulations of a model for 7-Fe2C>3 (maghemite) single particle of spherical shape 
are presented aiming at the elucidation the specific role played by the finite size and the surface 
on the anomalous magnetic behaviour observed in small particle systems at low temperature. The 
influence of the finite-size effects on the equilibrium properties of extensive magnitudes, field coolings 
and hysteresis loops is studied an compared to the results for periodic boundaries. It is shown 
that for the smallest sizes the thermal demagnetization of the surface completely dominates the 
magnetization while the behaviour of the core is similar to that of the periodic boundary case, 
independently of D. The change in shape of the hysteresis loops with D demonstrates that the 
reversal mode is strongly influenced by the presence of broken links and disorder at the surface. 

PACS Numbers: 05.10 Ln, 75.40 Cx, 75.40.Mg, 75.50 Gg, 75.50 Tf, 75.60 Ej 



7-Fe2C>3 (maghemite) and other magnetic oxides in the 
form of nanometric particles display anomalous magnetic 
properties at low temperatures. Experiments have shown 
that the hysteresis loops display high closure fields and do 
not saturateErB even at fields of the order of 50 T. Low 
magnetization as compared to bulk, shifted loops after 
field cooling and irreversibilities between the field cool- 
ing and zero fielcLcooling processes even at high fields 
are also observedoa. Moreover, the existence of aging 
phenomenaa in the time-dependence of the magnetiza- 
tion, indicates that there must be some kind of freezing 
leading to a complex hierarchy of energy levels. Whether 
these phenomena can be ascribed to intrinsic properties 
of the particle itself (spin-glass (SG) state of the surface 
which creates an exchange field on the core of the parti- 
cle) , or they are due to a,-coj.lective behaviour induced by 
interparticle interactionsBLj, has been the object of con- 
troversy in recent years and up to the moment there is no 
model giving a clear-cut explanation of this phenomenol- 
ogy, although simulation results for general small particle 
systemscl and, in particular, for maghemiteEl have been 
recently published. In order to contribute to elucidate 
this controversy we present the results of a Monte Carlo 
(MC) simulation of a single spherical particle which aims 
at clarifying what is the specific role of the finite size 
and surface on the magnetic properties of the particle, 
disregarding the interparticle interactions effects. 

Maghemite is a fcrrimagnctic spinel in which the mag- 
netic Fe 3+ ions with spin 5/2 are disposed in two sub- 
lattices with different oxigen coordination (8 tetrahedric 
(T) and 16 octahedric (O) sites per unit cell). In our 
model, the Fe 3+ magnetic ions are represented by Ising 
spins Sf = ±1 which allows to reproduce a case with 
strong uniaxial anisotropy while keeping computational 
efforts within reasonable limits. The spins interact via 
antifcrromagnetic exchange interactions with the near- 
est neighbours on both sublattices and with an external 
magnetic field H . In the simulation, we have used the 



reduced field h = [iHjkri in temperature units, where fi 
is the magnetic moment of the Fe 3+ ions. The values of 
the_nearest neighbour exchange constants for maghemite 
aret] J TT = -21 K, J 00 = -8.6 K, J TO = -28.1 K. We 
have used periodic boundary (PB) conditions to simulate 
the bulk properties and free boundaries (FB) for a spher- 
ically shaped particle with D unit cells in diameter when 
studying finite size effects. In the latter case, two differ- 
ent regions are distinguished: the surface formed by the 
outermost unit cells and an internal core. The size of the 
studied particles ranges from D — 3 to 14 corresponding 
to real particle diameters from 2.49 to 11.62 nm. 

We start by studying the effect of FB conditions and 
finite-size effects on the equilibrium properties in zero 
magnetic field. The simulations have been performed us- 
ing the standard Metropolis algorithm. Starting from a 
very high temperature (T — 200 K) and an initially disor- 
dered state with spins randomly oriented, the system was 
cooled down at a constant temperature step 8T = — 2 K. 
After discarding the first 1000 MC steps, the averages 
of the thermodynamic quantities were computed at each 
temperature during a number of MC steps ranging from 
10000 to 50000 depending on the system size. 

In Fig. [j], we compare the thermal dependence of the 
magnetization for spherical particles of different diam- 
eters D with the corresponding results for a system of 
size AT = 14 and PB (uppermost curve). A second order 
transition from paramagnetic to ferrimagnetic order sig- 
naled by a sharp peak at T C (D) in the susceptibility (see 
the Inset in Fig. [l]) is clearly observed. Finite-size effects 
on both the magnetization and the susceptibility are very 
important even for D's as large as 14 in the FB case, with 
T C (D) increasing as D increases and tending to the value 
for PB (which varies from 122 to 126 K when increasing 
N from 3 to 14) . The main feature in Fig. [l] is the reduc- 
tion of the magnetization M with respect to the PB case 
(dashed line) due to the lower coordination of the spins 
at the surface, which hinders perfect ferrimagnetic order 
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FIG. 1. Thermal dependence of the magnetization M. 
The results for particle diameters D = 3, 4, 5, 6, 8, 14 (from 
the lowermost curve in circles) and PB conditions N = 14 
(uppermost curve) are shown. Inset: Thermal dependence of 
the susceptibility for the same cases (the PB case is drawn in 
dashed lines). Munc is the ratio of the difference of O and T 
spins to the total number of spins. 

at finite temperatures. It is worthwhile to note that for 
all the studied diameters there is a temperature range, 
in which the demagnetization process of M is linear, this 
range being wider as the particle size decreases. In this 
linear regime the particle demagnetization becomes dom- 
inated by the surface effects, being the core and surface 
behaviours strongly correlated. Linear demagnetization 
is indicative of the effective 3D-2D dimensional reduction 
of the surface spips and has previously been observed in 
thin film systemst£l. 

Deeper insight on the magnetic ordering of this sys- 
tem can be gained by studying the thermal dependence 
of the equilibrium magnetization in a magnetic field. Sev- 
eral such curves are presented in Fig. ^|. They have been 
obtained by the same cooling procedure as described pre- 
viously in presence of different fields hpc- In this figure, 
the surface (continuous lines) and the core (dashed lines) 
contributions to the total magnetization have been dis- 
tinguished. For PB all the curves tend to the ferrimag- 
netic order value (i.e. M = 1/3). A maximum appears 
at high enough cooling fields hpc = 100 K which is due 
to the competition between the ferromagnetic alignment 
induced by the field and the spontaneous ferrimagnetic 
order (as the temperature is reduced the strength of the 
field is not enough as to reverse the spins into the field di- 
rection). However, for particles of finite size the curves at 
different fields do not converge to the ferrimagnetic value 
at low T, reaching higher values of the magnetization at 
T = the higher the cooling field (see lines with circles 
in Fig. ^|a). The total magnetization for small parti- 
cles is completely dominated by the surface contribution 
and this is the reason why the ferrimagnetic order is less 
perfect at these small sizes and the magnetic field can 
easily magnetize the system. However, the behaviour of 
the core is still very similar to that of the case with PB, 
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FIG. 2. Thermal dependence of M after cooling under 
Iifc- (a) corresponds to a spherical particle with D — 6. The 
results for two cooling fields Hfc — 20 K (lower curve) and 
Kfc = 100 K (upper curve) are shown. The contributions 
of the surface (thick lines) and the core (dashed lines) to the 
total magnetization (circles) have been plotted separatedly. 
(b) shows M for an N — 8 system and PB after cooling in 
Hfc = 20,40,80,100 K (from lower to uppermost curves). 
(See Fig. 1 for the definition of Munc)- 

although its contribution to M is small. At low fields, 
the surface is in a more disordered state than the core: 
its magnetization lies below M at temperatures for which 
the thermal energy dominates the Zeeman energy of the 
field. In contrast, a high field is able to magnetize the 
surface easier than the core due to the fact that the bro- 
ken links at the surface worsen the ferrimagnetic order 
while the core spins align towards the field direction in a 
more coherent way. 

In Fig. |p, we show the hysteresis loops for several 
particle diameters at 20 K. The loops have been com- 
puted by starting from a demagnetized state at h = 
(the first magnetization curve is not shown in this fig- 
ure). The results have been averaged over several (« 10) 
independent runs starting with different random num- 
ber seeds. First of all, let us note that the saturation 
field and high-field susceptibility increase as the particle 
size is reduced, since these quantities are mainly asso- 
ciated to the progressive alignment of the surface spins 
towards the field direction. As a consequence, the cycles 
of the smallest particles are similar to those observed in 
disordered systems, increasing their squareness (uniform 



2 




0.3 



0.0 



-0.3 



-0.6 







1 / / 

If / 

J 1 
J 1 


(b) 
i i 



-100 -50 50 
h(K) 



100 



FIG. 3. (a) Hysteresis loops for particles of diameters 
D = 3,4, 6, 8, 10 (from inner to outermost curves) and PB 
(long-dashed lines) at T — 20 K. (b) Surface (continuous line) 
and core (dashed line) contributions for a particle of diameter 
D = 3 at T = 20 K. 



rotation of M) with the size. In fact, by plotting sepa- 
ratedly the contributions of the core and surface to the 
total magnetization (see Fig. |l|b), we see that the loop of 
the core is almost perfectly squared independently of the 
particle size, indicating a coherent reversal of its magneti- 
zation, while the loop of the surface reveals a progressive 
reversal of M, which is a typical feature of a disordered 
or frustrated system. Nonetheless, the coercive field of 
the core is slightly higher but very similar to the one 
of the surface indicating that the reversal of the surface 
spins triggers the reversal of the core. We have shown 
that the existence of lower coordination at the surface of 
the particle hinders perfect ferrimagnetic order, increas- 
ing the magnetic disorder at the surface layer as the size 
of the particle decreases. However, the magnetic frus- 
tration associated to the competition between intra- and 
inter-sublattice antiferromagnetic interactions, finite-size 
and surface effects are not enough as to produce a SG 
layer contrary to the experimental observation of some 
authorsEm. This indicates that other ingredients (i.e., en- 
hanced anisotropy at the surface) should be included in 
the model to account for the SG layer. Our model shows 
that magnetic disorder at the surface simply facilitates 



the thermal demagnetization of the particle and also in- 
creases the magnetization at moderate fields, since sur- 
face disorder diminishes ferrimagntic correlations within 
the particle. 
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